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Abstract: 

We study eigenstates of chaotic billiards in the momentum representation and propose the radially 
integrated momentum distribution as useful measure to detect localization effects. For the momentum 
distribution, the radially integrated momentum distribution, and the angular integrated momentum 
distribution explicit formulae in terms of the normal derivative along the billiard boundary are derived. 
We present a detailed numerical study for the stadium and the cardioid billiard, which shows in several 
cases that the radially integrated momentum distribution is a good indicator of localized eigenstates, 
such as scars, or bouncing ball modes. We also find examples, where the localization is more strongly 
pronounced in position space than in momentum space, which we discuss in detail. Finally applications 
and generalizations are discussed. 
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1 Introduction 



In quantum mechanics the state of a system is given by a normalized vector in a Hilbert space, 
which in turn can be represented in different ways. Usually one chooses the representation 
which is most convenient for the problem at hand. In quantum chaos the main concern is to 
understand the semi classical limit, and the fingerprints which properties of the classical limit, 
like chaoticity or integrability, leave on the quantum system. 

A particulary convenient class of representations for the study of the quantum to classical 
correspondence is given by the Wigner function and its relatives, e.g. the Husimi density. The 
Wigner function of a quantum state is the representation which comes closest to a probability 
density on the phase space of the corresponding classical system [|l], || , and therefore should be 
a sensitive detector for classical fingerprints. One drawback of the Wigner function is that it 
is difficult to visualize, because it is a function on the 2n-dimensional phase space. Therefore 
one is often forced to study the projections on position or momentum space when one wants 
a visual representation of the state. Whereas the position representation is commonly used 
in quantum chaos, the momentum representation is not. Our aim here is to promote the 
momentum representation, by showing that it has some potential advantages. 
What we are especially looking for are fingerprints of the classical phase space structure, e.g. 
periodic orbits. So in the position representation one is searching for an enhanced probability 
density around an orbit, in order to detect e.g. so-called scarred eigenstates ||. But the 
detection of such a state from the probability density in position space can be rather ambiguous, 
there is no clear borderline which distinguishes scarred eigenstates from those which are not. 
In this case the probability density in momentum space can be helpful. E.g. assume we are 
studying an Euclidean billiard system, where the periodic orbits consist of segments of straight 
lines, therefore the momentum distribution of a state which is scarred by a periodic orbit should 
be concentrated on the directions of the line segments of that orbit. This criterion can be further 
simplified if one realizes that the momentum distribution of the n'th state is concentrated on 
the energy shell p 2 = E n , so that the main information is contained in the dependence on the 
direction of p. By integrating over the radial direction one gets a one-dimensional distribution, 
the angular distribution of the momentum, which contains all the essential information about 
the momentum distribution. This is the main advantage of the momentum representation, 
instead of working with a two-dimensional distribution as in position space, one can reduce the 
study to the one-dimensional angular distribution of the momentum. 

A related quantity was introduced in where the scattering approach is used to study ex- 
pectation values and scars on the Poincare section. The expectation values of a proposed "scar 
operator" correspond to a smoothed angular distribution of the momentum. Other approaches 
to study eigenfunctions consist in the use of the Bargmann representation and the stellar repre- 
sentation 0. For billiards these can be used to obtain representations on the Poincare section 
0, see also 0. 

The plan of the paper is the following. In the first section we use the boundary integral method 
to obtain expressions for the momentum representation of a given eigenfunction in terms of 
the normal derivative of this eigenfunction on the boundary. Furthermore we derive formulas 
for the angular and the radial distribution of the momentum density, respectively, expressed 
in terms of the normal derivative. In the second section we present a gallery of eigenstates 
displayed in position space and in momentum space, and compare the different representations. 
Finally applications and generalizations are discussed. 



2 Momentum distributions 



The boundary integral method is a common method for computing eigenvalues and eigenfunc- 
tions of two-dimensional billiards, see e.g. M, [|. The main point is that one can reduce the 
two-dimensional eigenvalue problem in such a billiard with the help of Green's formula and a 
Green function into an integral equation for the normal derivative of the eigenfunction on the 
boundary. Therefore one needs to solve only a one-dimensional problem, which is much more 
efficient numerically. 

2.1 Reduction to the boundary 

Let C I 2 be a domain in the Euclidean plane with piecewise smooth boundary dQ. The 
Dirichlet quantum billiard in Q is defined by 

-AMq) = PlMQ) for qeQ\dQ (1) 

ip n (q) = for q E d£l . (2) 

And we will use the trivial extension of ip n {q) to IR 2 by requiring ip„(q) = for q e IR 2 \ Q to 
M 2 . The transformation of this problem to a problem on the boundary rests on the formula of 
Green, 

J [f(q)Ag(q) - g(q)Af(q)] d 2 q = J [f(q(s))n(s)Vg(q(s)) - g(q(s))n(s)V f(q(s))] ds , (3) 

n an 

where / and g are arbitrary smooth functions, and n(s) denotes the outer normal on dQ, which 
is defined almost everywhere if the boundary is piecewise smooth. If one inserts for / a Green 
function, e.g. f(q) = \H^\p n \q — q'\), where Hq = Jo + i^o is the Hankel function of the first 
kind, which satisfies 

(A q + pl) 1 -Hl ) 1 \p n \q-q'\)=5(q-q>) , (4) 
and for g the eigenfunction x/) n , one gets for q el 2 \ dQ 

^Pn{q) = J ^H^(p n \q(s) - q\)u n (s) ds , (5) 
an 

where u n (s) := n(s)'Vip n (q(s)) denotes the normal derivative of ip n on the boundary expressed 
in terms of the arclength coordinate s. Eq. (^|) is the crucial relation by which one recovers the 
eigenfunction from its normal derivative on the boundary. Taking the normal derivative of (|5|) 
leads to an integral equation for u n on the boundary, which is then solved numerically. The 
choice of Hq for the Green function is necessary in order that this integral equation for u n 
does not lead to any spurious solutions, see e.g. ||. For the integral @ the contribution of Jq 
is irrelevant, because with (A + p 2 ) Jo(p n \q\) = the Green formula gives for q e Q \ dfl 

J Jo(Pn\q(s) - q\) u n (s) ds = , (6) 
an 



and so one has the simpler formula for the eigenfunction for q £ Q \ dQ 

i>n{q) = ~\f Y o(Pn\Q( s ) - Q\)u n {s) ds . (7) 



an 



Furthermore, one can express the L 2 -norm of ip n by aim 



^ J n(s)q(s)\u n (s)\ 2 ds=\\^ n \\ 2 . (8) 
n an 

This expression is very convenient as it allows to obtain a normalized normal derivative u n 
such that all other derived quantities of interest, such as expectation values, are correctly 
normalized. This relation is much faster to compute numerically than a normalization obtained 
from integrating |?/;(qr)| 2 over the billiard domain Q. 

Even in cases where the eigenf unctions are given by an expansion into some basis it can be 
useful to perform the computations starting from the normal derivative, as one can exploit the 
relation (j8|) to obtain normalized eigenfunctions. 

2.2 Momentum representation of eigenfunctions 

The momentum representation of an eigenfunction ip n is given by its Fourier transform 

Mp) :=^f e-™Mq) d2 Q = ^J e~ ipq Mq) d 2 g • (9) 

Using Green's formula (§) with f(q) = ip n (q) and g(q) = ■^ e ~ 1QP one gets 
! \ihJn)(-\r}) 2 p- illp - p- iqp (-r> 2 UJ 

2tt 



l - J [^(g)(-ip) 2 e-^-e^(-^)^( g )] d 2 g 

= hj [^(9(^))(-^(s)p)e- i9Wp -Wn(s)e- i ^] ds. 



(10) 



As ip n {q) = for q G dfl we get a representation of the Fourier transform as an integral over 
the normal derivative at the boundary 

MP) = I u n (s)e-^ p ds . (11) 



P 2 ~ Pi 2tt 

an 

This expression has an apparent singularity at p 2 = p 2 , but since the Fourier transform and its 
derivatives are bounded this singularity has to be cancelled by a zero of the integral. Solving 
fllTf) for the integral gives 

i- J u n (s)e-^ p ds = {p 2 - p 2 n ) Mp) , (12) 
an 

and taking the derivative of this equation with respect to r := |p| at the point r = p n one gets 

MP) = --^ J e~ ip ^pq(s)u n (s) ds , (13) 



an 



for |p| = p n . By repeating this procedure one can get expressions for the derivatives of ip n {p) 
at \p\ = p n as well. This can be used to interpolate near \p\ = p n for a numerical plot of ip n (p). 



2.3 Radially integrated momentum distribution 



In a billiard the classical flows at different energies p 2 are isomorphic, because the system scales 
with energy. At the quantum mechanical side the momentum distribution of the n-th state is 
concentrated around the energy shell p 2 = p 2 n . Because of the scaling property of the classical 
system the interesting phenomena, i.e. those which are due to the special system at hand, 
should occur only in the angular distribution of the momentum-probability. 
Therefore we want to study the angular distribution of the momentum defined by 

oo 

Info) / i>n(r,<p) r dr , (14) 



where we have introduced polar coordinates p = (r cos ip , r sin ip) . The advantage of this quan- 
tity, in contrast to the full momentum probability density or the probability distribution in 
position space, is that it depends only on one variable (p G [0, 2tt}. 
By inserting the expression fllT]) for the Fourier transform ip n into ( JH| ) one gets 



^(f) = 7^)2 j ( r 2_ p 2)2 J j e ^ Un{s)u n { S ') ds ds' dr , (15) 

n dUxdU 

where we have used the abbreviation 

a = a(<p,s,s'):=p(<p)(q(s)-q(s')) , (16) 

and p((p) = (cos cp, sin cp) is the unit vector in p direction. Using pj^p = — r 2 — p 2 anc ^ 
partial integration one obtains 

oo 

Art 2 I n (ip) = -2n 2 p 2 n |^ n (0)| 2 - ~ J -^-^ J j iae- [ra u n (s)u n (s') dsds' dr , (17) 

o n anxdn 

where it was used that J m u n (s) ds = — 2irp 2 l 'il) n (0) which follows from (p|). 
Here we have to discuss a problem which always occurs when inserting the expression (|TT| ) 
for ip n in an integral such as in eq. fll4l) . Since (|ll|) is the product of a factor which becomes 
singular at \p\ = p n and a factor which is zero there, we have an apparent problem if we want 
to interchange the order of integration in (|i~7D . Therefore we will choose a regularization for 
the factor 2 ]_ 2 . This can be done by adding to p n a small positive or negative imaginary 

P Pn 

part, p n ± ie, performing all computations and taking afterwards the limit e — 0. The result is 
independent of the regularization. In the following we will choose a symmetric regularization, 
i.e., replace (p 2 — p 2 ,)^ 1 by 

2 {p 2 - ( Pn + ie) 2 + p 2 - ( Pn - ie) 2 ) ' (18) 

which corresponds in the limit e — * to taking the principal values of all integrals over r = \p\. 
To avoid complicated formulas we will proceed formally without regularization, and interchange 
integrals whenever needed, but interpret all integrals over r as principal values. The reader can 
always check that the described regularization leads to the same results. 



Returning to the computation of I n (f) we interchange the order of integration in (|T7|) and get 

Iniv) = J J K n (ip,s,s')u n (s)u n (s') dsds' , (19) 



where the kernel is given by 



K n (tp,s,s') 



dQxdn 



dr 
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2p 



(20) 



a 



8tt 2 



sin(ra) 1 a 
7 dr - r-^r - i- — 



cos(ra ) 
r 2 -p 2 n 



dr 



and the integrals are interpreted as principal values. Note that the real part is an even function 
of a and the imaginary part is an odd function of a, thus K n has the property 



K n ((p,s,s') = K n (ip, s',s) 



(21) 



which is equivalent to I n {f) = In{f)- For Euclidean billiards the eigenfunctions can be choosen 
to be real. In this case u n (s)u n (s') is even under interchange of s and s', and therefore the 
integral over the imaginary part of K n , which is odd, vanishes. 
From [11]] one gets for the integrals (interpreted as principal values) 

oo 

1 



a. 



sm(ra ) 
r 2 -pi 



dr 



8tt 2 



(sin(|a|p n ) Gi(|a|p n ) - cos(|a|p n ) Si(|a|p n )) 



(22) 



a 



cos (r en J 

rn 



dr 



a 
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sin(|a|p r , 



(23) 



where Si and Ci denote the sine and cosine integrals. Note, that for billiards with Dirichlet 
boundary conditions p n > 0. With the abbreviation 

f(x) := sin(x) Ci(x) — cos(x) Si(x) , (24) 

we can write our final expression for I n (f) in terms of the normal derivative u n as 



In(<p) 



anxdn 



u n (s) u n (s') ds ds' 



+ i 

dnxdQ. 



a 



IGnp, 



sm(\a\p n ) u n (s) u n (s') dsds' . 



(25) 



Recall that a = a(<p,s,s') = p(ip)(q(s) — q(s')). If the normal derivative u n is real then the 
second term in (|2~5 ) vanishes, and if furthermore the eigenfunction ip n is odd with respect to 
some discrete symmetry of the billiard, see the appendix, then J Qn u n (s) ds = 0. So under these 
special conditions the expression ( [25] ) simplifies to 



Info) 



f(\a\p n )u n (s)u n (s') dsds' . 



(26) 



dflxdn 



2.4 Angular integrated momentum density 



Similarly to the radially integrated momentum density one can study the radial distribution of 
the momentum density by integrating over the angle (p. The radial momentum distribution is 
defined as 



2tt 



RJr) := r 



ipn{r, <p) 



dp 



(27) 



and if we insert the expression (|TTD for ip n one gets 

r 1 
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Rn(r) : 



Ur. 



n{s)Un{s') / e -Ms)-l(s')\rcos( V ) d(p dg dg > 



dflxdn 



(r 2 — p 2 ) 2 2ir 



(28) 



anxdn 



Again one has an apparent singularity at r = p n due to the pre-factor l/(r 2 — p^) 2 - In the 
same way as in the previous section one obtains by differentiating (r 2 — p^) 2 R n (r) / r sufficiently 
often 



-(- 

8p n \dr 



R 'n(.Pr, 



24p n V dry 



d \' 



/ 2 2 \ 2 - 1 L n 
i r -Pn) 



Rn{r) 



RJr) 



r=p n 



r=p n 



(29) 
(30) 



and higher derivatives of R n at r = p n can be obtained in the same way. Using eq. (|28|) one 

gets 



Rn{Pn) 



1 



167Tp r , 



u n (s) u n (s')\q{s) - q{s')\ 2 



dflxdn 



[J 2 (\q(s) - q(s')\r) - J (\q(s) - q(s')\r)} dsds' (31) 



R'niPr, 



A8np r , 



u n {s)u n {s )\q{s) - q{s 



dflxdn 



- [3U\q(s) - q(s')\r) - J 3 (\q(s) - q(s')\r)} dsds' , (32) 



and these expressions can be used to interpolate in the region around r = p n . 



3 Gallery of eigenfunctions in momentum space 



In the following we present a number of examples of eigenfunctions in momentum representa- 
tion for two chaotic billiards and compare them with the position representation. A series of 
eigenfunctions of the cosine billiard in momentum representation can be found in |12 |. The 



first system we study is the stadium billiard, which is proven to be strongly chaotic, i.e. it is 
ergodic, mixing and a .fT-system fi~3| , |i"4 |. The height of the desymmetrized billiard is chosen 
to be 1, and a denotes the length of the upper horizontal line, for which we have a = 1.8 in the 
following. To compare the structures in the eigenfunctions with the classical orbits, we use the 



symbolic dynamics proposed in [15], which is proven in [16 



The second system is the cardioid billiard, which is the limiting case of a family of billiards 
introduced in ||17|| . The cardioid billiard is proven to be ergodic, mixing, a i^-system and a 
Bernoulli system [18|, O, |20[|. The eigenvalues of the cardioid billiard have been provided by 



see e.g. [|2|, |3 



proposed in |2J, |25 



Prosen and Robnik |^1[ and were calculated by means of the conformal mapping technique, 
To describe the periodic orbits in the cardioid we use the symbolic dynamics 
:e ]Tj], [26j for proofs. 

As these two billiards are ergodic, the quantum ergodicity theorem |27], |28], ^9], [K], |31| applies, 
see e.g. 



for introductions. It implies that for "almost all" eigenfunctions their Wigner 



32 



functions W n (p, q) become equidistributed on the energy shell in the semiclassical limit, i.e. 



W ni (p,q) :-- 



3 i<7 P 



Q- 



^n, q + 



2„/ 



d z q 



(33) 



where vol(Se) = ff R 2 xn $(p 2 — E) d 2 pd 2 g is the volume of the energy shell. "Almost all" 
precisely means that eq. (^) holds for a subsequence {rij} C N of density one, that is one has 



#{n ? - E n . < E} 
lim ,\ 3 " J — — - = 1 

e^oo #{n | E n < E\ 



(34) 



For the eigenstates in position and momentum reperesentation the quantum ergodicity theorem 
implies that 



IVv(<z)p 



vol (ft) 



as n, 



oo 



and 



1 



oo 



I "0n, (p) 1 2 -> - $ (P 2 ~ E n ) as rij 
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in the weak sense. For the radially integrated momentum distribution it follows that 

J n f {<p) 



1 

2t7 



as rij 



oo 



(35) 



(36) 



(37) 



in the weak sense. For the angular integrated momentum density one has in general, indepen- 
dent of the ergodic properties of the classical system, 



Rr, 



2r5(r 2 - E n ) as rtj 



(38) 



Of special interest are of course the subsequences of exceptional eigenfunctions (if they exist) 
which are according to the quantum ergodicity theorem of density zero for ergodic systems. A 
drastic example of such eigenfunctions are the so-called bouncing ball modes in the stadium 
(which also occur in other billiards with two parallel walls). These eigenfunctions localize on the 
orbits which bounce up and down between the two parallel walls with perpendicular reflections 



at the boundary. For the stadium billiard the counting function for these modes increases 
asymptotically as E 3 ^ 4 |35|, |36| . 

Another important class of exceptional eigenfunctions are so-called "scarred" eigenfunctions 
H, showing localization along unstable periodic orbits, which have been first observed in the 
stadium billiard. Also in the cardioid billiard scarred eigenfunctions have been observed [p4 . 37[ . 



As has been mentioned in the introduction, for an eigenfunction which is scarred by a periodic 
orbit, Ini'fi) should have prominent peaks at the angles corresponding to the directions of the 
orbit, whose intensity is expected to be proportional to the length of the corresponding orbit 
segment. 

Let us start the discussion by showing plots of some low lying eigenfunctions of the stadium 
billiard with odd-odd symmetry. Figs. |l]-|2| show for the stadium billiard three-dimensional plots 
of \ip n (q)\ 2 , |^n(p)| 2 and the corresponding grey-scale plots. Furthermore I n (ip) and R n {r) are 
shown. In the case of I n ((f) we only plot the interval ip G [0,7r/2], as the other directions are 
obtained from symmetry, e.g. I n (f) = I n {^ ~ V 9 ) f° r V 9 £ [ 7r /2, 7r]- The dashed line in the plot 
of I n ((p) corresponds to the mean value l/(27r), see eq. d37|) . The state n = 24 displayed in 
fig. [TJ does not show any prominent localization in position space, whereas the state n = 26, 
fig. 0, is an example of a low lying bouncing ball mode. For both eigenstates the momentum 
distribution |^ n (p)| 2 is mainly concentrated around the energy shell E n = p 2 , which is indicated 
by the inner full circle in figs. This is also nicely seen in the plots of ^24(7") and Rieir) 
where the radius of the energy shell is marked by a rhombus. Whereas |V'24(p)| 2 shows peaks 
at several places, \^2&{p)\ 2 only has prominent peaks around the p y direction. The pictures 
shown in fig. are precisely those expected for a bouncing ball mode |3(J: in position space 



we have localization on the rectangular part of the billiard, and in momentum space there is 
localization in the direction of the bouncing ball orbits. That there are four major peaks of 
I ^26 (p) 1 2 visible in fig. |2] is due to the symmetry of the billiard, which implies that \ip n (p)\ 2 is 
reflection symmetric with respect to the axes p x = and p y = 0, on which |-?/> n (p)| 2 vanishes 
for the stadium billiard with everywhere Dirichlet boundary conditions. According to RH] the 
peaks should be in the directions 

i^g and vr ± (fl h t with := arctan(a//A;) , (39) 

where k denotes the number of modes in x-direction and I denotes the number of modes in 
^/-direction and a = 1.8 is the billiard parameter. If the ratio l/k increases, then the value 
of ip^ increases, and in the limit l/k — > oo we have (p^\ — » n/2. For the bouncing ball mode 
shown in fig. ^| there is good agreement with the direction ip^s = arctan(4a), which is marked 
by a full triangle in the plot of hei^p)- 

Fig. |3]a) shows an example of an eigenfunction having a quite regular and uniform pattern in 

position space. However, the plots of |^263(p)| 2 an d h&zi.v) reveal three major momentum 

directions. Although we did not succeed to find an orbit corresponding to these directions, 

there is some indication that this orbit has a quite large geometric length. 

The second eigenfunction shown in fig. |3] is an example of a higher lying bouncing ball mode 

with 14 modes in x-direction and 30 modes in y-direction. For this bouncing ball mode the 

distance of the maximum to these directions is larger than for the one displayed in fig. 0. This 

is well accounted for by the formula ( |3"9"D as the location of the maximum of /45s ((p) agrees well 

with v?i4 3o = arctan(30/14 a), which is marked by a full triangle in the plot of /45s(<^). 

Fig. [| shows two examples of eigenfunctions showing scarred structures in position space. 

1^1771 (?) | 2 shows an enhanced probability along the two orbits with symbolic codes 0123 for the 



orbit of diamond shape, and 1155 for the orbit of rectangular shape (shown as dashed lines in 
the plot of |^i77i(q)| 2 )- These two structures are less visible in the momentum representations. 
For the contribution of the 0123 orbit to Iim{ip) we find that it is spread out near to <p = tt/8, 
in agreement with the structure in position space, which is also not aligned precisely along the 
0123 orbit. The higher intensity near to ip = corresponds to localization along the 1155 orbit. 
Although \ipnn{q)\ 2 shows stronger localization near to the circular boundary which seems to 
be associated with the (p = ±rr/2 directions, there is no prominent peak in that direction. How- 
ever, there is some enhancement of the total probability to find the particle with momentum 
directions 3tt/8 < <p < rr/2 which is seen more clearly in the plot of |^mi(p)| 2 than in the plot 
of 

Imiitp)- The rapid oscillations of Iim(<p) near to ip = rr/2 are not a numerical artefact; the 
structures are also visible in the plot of \ip\m {p)\ 2 i we wm return to this point at the end of this 
section. The eigenfunction \ipm6(q)\ 2 displayed in fig. |b) shows localization along the 1133 
orbit. From this one would expect two major directions in momentum space. However, this is 
not seen in the pictures of |?/>i776(p)| 2 and Imei'p)', one observes several significant momentum 
directions, and at first sight no simple relation to the plot of 1^1776 (^) | 2 seems to exists. But 
there are three peaks in ii776(y?) close to one of the directions of the triangular orbit, so there 
is an enhanced probability around the orbit in momentum space. This seems to be a common 
phenomenon, when interpreting the pictures of I n {<p), one is sometimes mislead by the height 
of some peaks, the relevant information is the area under the curve around this peak, which 
can be quite small due to the narrowness of some of the peaks. 

Fig. |5] shows two eigenfunctions with clear localization in position space, this time for the 
stadium billiard with even-even symmetry. The first is localized along the unstable 01 orbit 
running along the symmetry-axis. This localization is also clearly seen in momentum space, 
where \^>5oo{p)\ 2 shows strong enhancement in the directions tp — 0, 7r, corresponding to a strong 
peak of hoo(f) at <p = 0. In case of the eigenfunction displayed in fig. ||b) there is also just one 
distinguished momentum direction. The orbit 111333 has two momentum directions, where 
the one with smaller angle is expected to dominate due to the longer line segments. This is 
also seen in the plot of /1273 (</?), where the momentum direction corresponding to the shorter 
segments of the orbit do not acquire a high probability. 

Fig. H shows two further examples of localized eigenfunctions. The first has a nice pattern in 
position space, which appear to be associated to the 10203020 orbit (only one orbit of the two 
symmetric partners is shown in fig. ^a). This orbit has just one momentum direction, however 
in momentum space we observe five prominent momentum directions. One could associate the 
three leftmost peaks in Iuoq(p) to this orbit, the additional peaks might be caused by some fine 
structure, which is hardly visible in the plot of | ^1409(9) | 2 - The eigenfunction with n = 1993 
displayed in fig. also shows clear localization in position space along the 110552 orbit. 
In momentum space there is one significant momentum direction at <p m 37r/8, however this 
direction does not correspond to the pattern seen in position space, from which one would expect 
an important contribution to | ^1993(2?) | 2 and ii99 3 (<^) in the direction of <p = 0.71 .... Indeed 
near to this direction there is a range of angles for which /1993 (y) has a number of smaller peaks, 
which together lead to an enhanced probability in this interval compared to other intervals. 
This example again illustrates that not just the height of the peaks is of importance in the 
interpretation of the pictures for I n (<p), but the overall accumulated probability corresponding 
to some interval of momentum directions. 



In [34] the rate of quantum ergodicity has been studied both in position and in momen- 



tum space. For the stadium billiard it was expected that if the bouncing ball modes are 
the dominating (in the sense of having the strongest increase in the counting function) sub- 



n = 24, odd-odd symmetry 




Figure 1: Three-dimensional plots of 1^24(9) \ 2 , 1^24 (p)| 2 , their corresponding grey- 
scale pictures and the plot of the radially integrated momentum distribution ^a^) 
and the angular integrated momentum distribution R2i{f)- The momentum distri- 
bution \ip2i{p)\ 2 is concentrated around the energy shell, which is indicated as the 
inner circle. This is also clearly visible in the plot of i?24( r )- The angular distri- 
bution I2a{<^) does not show any significantly preferred directions and the plot of 
^24(9) | 2 also does not show any prominent patterns. 



n = 26, odd-odd symmetry 




71/4 cp 71/2 5 10 15 20 r 25 

Figure 2: Example of a low-lying bouncing ball mode, for which \ip26{q)\ 2 is localized 
on the rectangular part of the billiard. The plot of the momentum distribution 
|-?/> 26 (p)| 2 shows a strong localization in the momentum directions (p x ~ 0, p y ~ ±p n ) 
of the bouncing ball orbits. This is also clearly seen in the plot of -^(v 9 )) which is 
concentrated near to (p = n/2. The corresponding direction yj^s = arctan(4a) is 
marked by a triangle. In the plot of R26{r) one observes some additional oscillations 
in comparison to R2a{t) in fig. [I]. 



a) n — 263, odd-odd symmetry 



b) n — 455, odd-odd symmetry 






Figure 3: In a) an example of an eigenfunction is shown which appears to be com- 
pletely non-localized in position space. However, in momentum space 1^263 (p) | 2 
and /263(y) show clear localization in three major momentum directions. The sec- 
ond eigenfunction is an example of a higher lying bouncing ball mode with 14 modes 
in x direction and 30 modes in y direction. Consequently the corresponding mo- 
mentum distributions show localization near to the ip = n/2 direction. In the plot 
of l45s(ip) the triangle marks the direction ^14,30 — arctan(30/14 a). 



a) n — 1771, odd-odd symmetry 



b) n — 1776, odd-odd symmetry 




Figure 4: In a) the eigenfunction shows localization along the 1155 orbit of rect- 
angular shape and also the 0123 orbit in the stadium. In the plot of Inn((p) the 
corresponding momentum directions are marked by full triangles (1155 orbit) and 
open triangles (0123 orbit). The eigenfunction in b) shows scarring along the 1144 
orbit, which is clearly seen in the plot of |^i776(g)| 2 . However, in this case the scar- 
ring is not that clearly reflected in the pictures of |^i776(p)| 2 and I 1776 ({p), see the 
text for a discussion. 



a) n — 500, even-even symmetry 



b) n — 1273, even-even symmetry 
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Figure 5: For n = 500 one has a nice example of an eigenfunction localized around 
the shortest unstable orbit, which runs along the symmetry axis. This is also seen in 
the plots of IV'soo (j*) | 2 an d hooiv)^ which show a strong enhancement in the ip — 0, 7r 
directions. The eigenfunction in b) shows localization around the 111333 orbit. 



a) n — 1409, even-even symmetry 



b) n = 1993, even-even symmetry 
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Figure 6: The eigenfunction with index n = 1409 of even-even symmetry is localized 
along the 10203020 orbit (only one of the two symmetry-related partners is shown). 
Surprisingly the plots of | ^1409 (p)| 2 and ii409(y?) show five peaks instead of just one, 
see the text for further discussion. The eigenfunction in b) is localized in position 
space along the 1105552 orbit. This orbit has two different momentum directions (in 
the intervall ip G [0, 7r/2]), still there is an additional high peak visible for Iiggs^). 
This illustrates that for I n ((p) not just the height of the peak is relevant, but the 
total area below a peak, which accounts for a high probability in the corresponding 
direction. 



a) n — 1652, odd-odd symmetry 



b) n — 1874, odd-odd symmetry 




Figure 7: Two examples of eigenfunctions which are localized on the rectangular part 
of the stadium billiard, which are not bouncing ball-modes. These eigenfunctions 
may be understood in terms of the plotted orbits, although in momentum space 
there are additional prominent momentum directions. This might indicate that 
both orbits contribute to each eigenfunction, see the text for further discussion. 



sequence of localized eigenf unctions, that then the rate of quantum ergodicity should obey 
nJe] ^2e u <e Aipn) — <j(A)\ ~ E 1 / 4 . However, it turned out that the rate is 



much slower |34|| , at least in the considered energy range covering the first 6 000 eigenfunctions. 
This was attributed to a considerable number of eigenfunctions, which show localization in posi- 
tion space in the rectangular part of the billiard, without being bouncing ball modes, see fig. 11 
of p3 |. Two examples are shown in fig. [7|. One sees that they are also localized in momentum 
space near the direction <p = n/2, like the bouncing ball modes, but have more peaks. These 
states seem to correspond to the set of periodic orbits which are bouncing m times between the 
parallel walls before they are reflected into themselves in the circular part. Two of them, the 
202023202024 and the 2020202320202024 orbit, are shown in the density plots of |t/>i 6 52(<?)| 2 
and 1^1872 (q) | 2 - Since the lengths of these orbits are close to being rationally dependent one 
can speculate, that the naive Bohr-Sommerfeld quantization leads to two eigenvalue-sequences 
which have some very close pairs of eigenvalues, for which the corresponding eigenfunction 
should be scarred by both orbits. The plots of /i652(v) and /i872 (y) indicate that this happens, 
because both have strong peaks at the same directions corresponding to the two orbits. 
The cardioid billiard possesses no parallel walls and therefore no bouncing ball modes; in this 
sense it is more generic than the stadium billiard. Figs. [8|-[^ show for the cardioid billiard three- 
dimensional plots of \i/j n (q) | 2 > |"0ri(p)| 2 and the corresponding grey-scale plots. Furthermore 
I n {f) and R n {f) are shown. The eigenstate of odd symmetry with n = 24 shows localization 
along the unstable orbit AB, wich is shown in the plot too. This localization is also seen in the 
plots of 1^24 (i^) | 2 and ^(v 9 )- I n the plot of Iu{(p) the corresponding direction of the orbit is 
marked by a full triangle at <p = rr/2. The plot of R2a{t) shows the expected localization (plus 
some oscillations) around the energy shell r = \p n \, indicated by the rhombus. In the plots 
of the Fourier transforms the energy shell is marked by the full inner circle. The eigenstate 
displayed in fig. || does not show any particular localization both in position and in momentum 
space. 

In fig. [10] the first example shows clear localization in position space along the orbit with 
code AAABBB . In momentum space the vertical direction (ip = vr/2) shows an enhanced 
probability, whereas the second direction, ip = 0.62 . . . , is not as prominent as one might 
expect from the picture of the eigenfunction in position space. The second eigenfunction shown 
in fig. [10] is an example of a non-localized eigenfunction, which is also nicely seen in the plots 
of the momentum distributions, which do not show any preferred momentum direction. 
Fig. |TT] shows two examples of eigenfunctions localized along the AB orbit, which is clearly seen 
in the plots of | ^isi7 (<?) 1 2 and | "02605 (q) | 2 . The latter eigenfunction is the one which was found 
in p8| to have the largest maximum norm ||^ n ||oo among the first 6000 eigenfunctions of odd 



symmetry. Also the corresponding momentum distributions reveal that the direction (p = it/2 
stands out. However, for the eigenfunction with index n = 2605, there are more oscillations 
visible in the plot of Inif), but still with a high probability near ip = tt/2. 
The first eigenfunction displayed in fig. O shows localization along the cusp orbit CAAC, and 
thus is an example of an diffractive scar. The enhanced momentum directions seen in Ii 908 ((p) 
are not precisely at the places of the directions of the orbit CAAC, but of the same quality 
of the agreement of the structure visible for l^igos (^) | 2 with the shape of the orbit. For the 
second example shown in fig. |T_2] one observes that | "01538 (9) 1 2 is clearly localized along the 



billiard boundary. Due to the non-convexity of the billiard there are no whispering gallery 
orbits as in the stadium billard. However, as found in in the cardioid there exist families of 
periodic orbits which accumulate in length. A candidate responsible for the structure visible in 
|"0 1538 (qr)| 2 is the family of orbits with code A n BABB. This type of orbit also accounts for the 



n = 24, odd symmetry 




Figure 8: Three-dimensional plots of 1^24(9) | 2 ; 1^24 (p)| 2 , their corresponding grey- 
scale pictures and the plot of the radially integrated momentum distribution ^a^) 
and the angular integrated momentum distribution i?24( r )- The momentum dis- 
tribution |V'24(p)| 2 is concentrated around the energy shell, which is indicated by 
the inner circle. This is also clearly visible in the plot of R2i{r)- Furthermore this 
state is to some extent localized along the AB orbit, leading to an enhancement of 
1^24 (p)| 2 near to <p = tt/2, 3ir/2, also seen in the plot of I^djf) near to the momentum 
direction tp — tt/2 (marked by a triangle). 



n = 27, odd symmetry 
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Figure 9: Same as in the previous figure but for n = 27. In this case there is no 
prominent localization neither in position nor in momentum space. 



a) n — 1970, odd symmetry 



b) n — 1277, odd symmetry 




Figure 10: The eigenfunction in a) shows localization along the AAABBB orbit. 
There is some enhancement in the corresponding momentum directions, in partic- 
ular near tp = 7r/2, and less significantly for the second direction near tt/4. The 
eigenfunction in b) is an example which appears to be quite delocalized both in 
position and in momentum space. The pictures look like those expected, according 
to the quantum ergodicity theorem, for a typical eigenfunction . 



a) n — 1817, odd symmetry 



b) n — 2605, odd symmetry 




Figure 11: Two examples of eigenfunctions which are localized along the AB orbit. 
In a) the localization appears to be more pronounced than in b). In particular the 
enhancement near to tp = n/2 is much less pronounced for ^mbiv) i n comparison 
to hsniv)- 



a) n — 1908, odd symmetry 



b) n — 1538, odd symmetry 




Figure 12: The eigenfunction displayed in a) shows localization along the cusp 
orbit CAAC, and thus is an example of a diffractive scar. The association of the 
corresponding momentum directions seems to be ambiguous. The second example 
shows an eigenfunction, which is mainly localized along the boundary. The family of 
orbits with code A n BABB appears to describe the observed pattern quite well. In 
the plot of Iisz%{(p) only the momentum directions of the part of the orbit running 
inside the billiard are marked by triangles. 



a) n — 1252, odd symmetry 



b) n — 1915, odd symmetry 





Figure 13: In a) a nice example of an eigenfunction showing scarring along the 
AABB orbit is shown. Near to the symmetry axis the eigenfunction density is 
spread out. In momentum space there is a clear enhancement, which is close to 
(although not as close as one might expect) the momentum direction of the orbit. 
The eigenfunction displayed in b) is localized along the triangular orbit AAB. Again 
there are some deviations of the pattern of the eigenfunction near the symmetry 
axis. Also, the momentum distributions do not show a clear enhancement in the 
momentum directions of the orbit. 



additional structures visible in the interior of the billiard. As for an orbit running along the 
boundary all momentum directions occur, one expects to see only the two further directions of 
the part of the orbit inside the billiard. Indeed, there is an enhanced probability for these two 
directions, which are marked by triangles in the plot of /153s (</?)■ One should remark, that there 
is some additional structure visible in | < 01538(q) | 2 which could correspond to the AAABBB 
orbit (see fig. |T0|) . 

We conclude our survey of eigenf unctions in momentum space with two further examples of 
localized eigenfunctions, see fig. |13|. The first eigenfunction is localized along the A ABB orbit, 
which is also seen in the corresponding momentum distributions. Near to the symmetry axis 
|"0i252(<?)| 2 the enhanced region is much larger, presumably due to interference effects. The 
eigenfunction with n = 1915 of odd symmetry is localized along the triangular orbit AAB. 
Near the symmetry line the structures of the eigenfunction deviate from the shape of the orbit, 
which could be caused by some contribution of the CAAC orbit, or again by interference 
effects. Both directions of the AAB orbit are also seen well in the pictures of the momentum 
distribution. 

For the cardioid billiard one generally observes a finer scale on wich I n fluctuates, compared 
to the stadium billiard. This is an effect of the stadium billiard, the long parallel segments of 
the boundary are reflected in the momentum distributions as a stretching of the structures in 
the p y direction. Therefore the radial integration acts as a smoothing except near <p = where 
accordingly I n (<p) fluctuates on a much finer scale. 



4 Summary, Applications and Discussion 

In this paper we have proposed the representation of eigenfunctions in momentum space as 
an important and useful complementary picture to the commonly used position representa- 
tion. In particular we have introduced the radially integrated momentum distribution I n {'~p)i 
which captures all essential features of the momentum distribution \ip n (p)\ 2 - As I n (^>) is just a 
one-dimensional function, it provides the information on the momentum distribution in a very 
condensed form. Explicit formulae for ip n (p) and I n {<p), and also the angular integrated mo- 
mentum distribution R n (r), are given in terms of the normal derivative u n (s) along the billiard 
boundary, which allows for an efficient numerical computation of these quantities. 
For the stadium and the cardioid billiard several examples of \ip n (q)\ 2 , \ipn(p)\ 2 , I n (^p), and 
R n {r) are given. For the stadium billiard in particular the bouncing ball modes lead to signifi- 
cant peaks in the radially integrated momentum densities I n (tp). Also eigenfunctions showing 
scarring in position space lead to scarred states in momentum space. In addition to eigenstates 
showing localization around unstable periodic orbits, we have found for the cardioid billiard ex- 
amples of eigenfunctions which are localized along cusp orbits, which run into the singularity of 
the cardioid billiard. Such diffractive scars are surprisingly clearly visible, despite the fact that 
the contribution of the corresponding orbit to the density of states is of lower order in h than 
the periodic orbits [0]. Our numerical studies also reveal that there are eigenstates, for which 



the correspondence between localization in position space and localization in momentum space 
is not as clear as one would expect. We see two main reasons for this: firstly, the fluctuations 
appear to be larger in momentum space than in position space (this is also seen in the results 
for the rate of quantum ergodicity in momentum space, which appears to be slower than in 
position space ]34|]). Secondly, scars which are weak both in position and momentum space are 



better visible in position space. This is because the usually symmetric pattern of the scar can 



be detected by the eye more easily against the fluctuating background in the two-dimensional 
plot of |^„(<?)| 2 , than in fc(p)| 2 or I n ((p). 

An application of the radially integrated momentum distribution I n (<p) could be to use it to 
construct a quantity which detects localized eigenf unctions. A suitable definition might be 
jscar ._ _L- f^^g I n ((p) dip, where U is the geometric length of the orbit segment corre- 
sponding to the direction (pi and 5 is some appropriate width, which might be chosen to be 
energy dependent. In this context it may be also useful to employ the results of the semiclas- 
sical studies of eigenfunctions |39|, [5], [41] to obtain expressions for I n ((p) in terms of periodic 



orbits. A number of different scar measures have been defined and studied, like for example the 
integral of the Wigner function over a tube around the orbit in phase space [41], the integral 



over a tube in position space [42 



sec 



44 for a 



or the use of quantum Poincare sections |43| . 
recent review and references therein. 

The use of momentum distributions as representations of eigenstates may also be useful for other 
systems, like systems with potential etc. The definition of the radially integrated momentum 
density may easily be generalized to other scaling systems. The usefulness of I n ((p) is due to 
the fact that the classical motion occurs on straight lines in Euclidean billiards. If this does 
not hold, the above simple interpretation of I n (<p) is lost. 

We believe that the use of momentum distributions, in particular the radially integrated mo- 
mentum distribution is a convenient and useful representation, providing additional information 
to the commonly used position space representation. 
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Appendix: Symmetry reduction 

When the billiard has some symmetries they can be used to reduce the domain of integration 
over the boundary to a smaller part in integrals as (|2lf). We will discuss as examples the 
stadium billiard and the cardioid billiard. To simplify the notation we use q = {x,y). 
The stadium billiard has two symmetries, it is symmetric with respect to reflection at the x-axis 
and at the ?/-axis. Therefore the eigenfunctions can be choosen to be even or odd with respect 
to each of these symmetries. So we will treat the case of a function ip(x, y) which satisfies 

ij)(-x,y) = (-1)S(^2/) (40) 
t/>(x,-y) = (-1)^(^,2/) , (41) 

with l x ,l y E {0, 1}. Thus ip either satisfies Dirichlet or Neumann boundary conditions on the 
symmetry axes. The symmetry of the wavefunction implies now certain symmetries of the 
normal derivative u(s) = n(s)(Vip)(q(s)). Denoting by L the total length of the boundary of 
the full system, we can express u(s) for s E [L/4,L/2], s £ [L/2, 3L/4] and s E [3L/4, L] in 



terms of u(s) for s G [0, L/4]. Here we chose s = to be the intersection of the horizontal 
symmetry axis with the right semicircle, and s is oriented counterclockwise. One gets for 
8 E [0, L/4] 

u(L/2 - s) = (-l) l *u(s) (42) 
u{L/2 + s) = {-l) l * +i »u{s) (43) 
u(L -s) = (-l) l *u(s) . (44) 

These relations can be used to reduce an integral over the boundary f(s)u(s) ds to an 
integral over one quarter of the boundary. For instance with ([E|) one gets 

L/2 o 

u(s)f(s)ds= J u(L/2 + s)f(L/2 + s) ds 

-V4 (45) 

L/4 L/4 

u(L/2 - s)f(L/2 -s)ds= [ {-l) l *u{s)f{L/2 - s) ds , 



and applying ( ^3|) and d44|) to the other parts of the integral leads to 

L L/4 



u(s)f(s)ds= j u(s)[f(s) + (-l) l *f(L/2-s) 



+ (-l) l * +l »f(L/2 + s) + (-l) h f(L - s)} ds 



(46) 



which is the desired reduction of the integral. We have to treat as well double integrals over 



the boundary of the form u(s) u(s') f(s, s') ds ds'. If we use ( fi6D we get an integral with 

16 terms, but for special f(s, s') this can be reduced further. If / is of the form 

f(s,s')=g(\ q (s)-q(s')\) (47) 

which is the case in our applications, then we have 

f(s, s') = /(L/2 - a, L/2 - s') = /(L/2 + s, L/2 + s') =f(L-s,L- s 1 ) , (48) 

because of the symmetry of the billiard the distance between the points on the boundary is the 
same for all four pairs of arguments of /. Applying ( |l~8f) and ([46D to a double integral over a 
function of the type ( f4?D leads to 

L L L/4 L/4 

u{s) u(s') f(s, s') ds ds' = [ [ u{s) u(s')f(s, s') ds ds' , (49) 



with 




oo oo 



f(s,s') =4[f(s,s') + (-l) l *f(s,L-s') 

+ (-l) l ^f(s, L/2 + s') + (-!)<* f(s, L/2 - s')] . 



The cardioid billiard has only one reflection symmetry at the x-axis, so the eigenfunctions can 
be chosen to be even or odd, 

il;{x,-y) = (-l) l »i)(x,y) (51) 
with l y e {0, 1}. For the normal derivative this leads to 

U {L - s) = (-l)St(s) , (52) 
and so we get for an integral over the boundary 

L L/2 

J u(s)f(s) ds= J u(s) [f (a) + (-l) l »f(L - s)] ds . (53) 



For a double integral over a function of the type (f|7]) one obtains 

L L L/2 L/2 

J Ju{s)u{s')f{s,s')dsds' = J J u(s)u(s')2[f(s,s') + (-l) l «f(s,L-s')} dsds' . (54) 
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